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Abstract 

We have studied a modified Yang-Mills-Higgs system coupled to Einstein gravity. The modification of 
the Einstein-Hilbert action involves a direct coupling of the Higgs field to the scalar curvature. In this 
£f*l • modified system we are able to write a Bogomol'nyi type condition in curved space and demonstrate that 

' the positive static energy functional is bounded from below. We then investigate non-Abelian spherically 

symmetric static solutions in a similar fashion to the 't Hooft-Polyakov monopole. After reviewing previ- 
ously studied monopole solutions of this type, we extend the formalism to included electric charge and we 
present dyon solutions. 

(N ' 
O 

1 Introduction 

In a classic paper Q| , Dirac proposed the possible existence of a magnetic monopole, the analogue of an isolated 
electrically charged particle. Motivated principally to restore the symmetries between electric and magnetic 
forces, Dirac found that the existence of a monopole provided a natural explanation for the quantization 
of electric charge. A description of such a monopole consistent with quantum mechanics would lead to the 
■ famous Dirac charge quantization condition, 

eg = -jnhc, 

where e 2 /Tic is the fine structure constant, g is the monopole charge, and n an integer. Dirac's theory 
required a U(l) valued gauge potential which was singular along a line (Dirac string) originating from the 
monopole and extending to infinity. Later Dirac's theory was reformulated by Wu and Yang Q within the 
framework of fiber bundles. The singular line is avoided at the expense of introducing coordinate patches on 
a sphere surrounding the monopole. However, the transition functions between the coordinate patches that 
are elements of the U(l) gauge group are singular. Thus, there is a complete equivalence between the two 
descriptions (^] and the monopole emerges on a sound footing like any other particle in nature. Consequently, 
during the past decades, extensive effort has gone into experimental search for monopoles, but unfortunately 
has as yet had no success. 

In spite of the lack of experimental success, the monopole continues to thrive in the theoretical laboratory. 
With the pioneering work of 't Hooft and Polyakov [Q , the monopole was reinvented in a new form as a finite 
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energy, particle-like soliton in non-Abelian gauge theories with spontaneous symmetry breaking. Moreover, 
such objects are generic in any spontaneously broken non-Abelian gauge theory which has an unbroken U(l) 
gauge symmetry. Such monopoles are expected to be produced in abundance in phase transitions of grand 
unified theories, which has implications for early universe cosmology. New searches for such relic monopoles 
are under way |||(|. 

More recently, a great deal of activity has centered around monopoles in curved space-time in order to 
study the effects of gravity. New insights have emerged from a study of a coupled Einstein- Yang-Mills-Higgs 
[EYMH] system with solutions describing black holes with magnetic charge, black holes within magnetic 
monopoles and magnetic monopoles within black holes fljO]. Consequently, gravity cannot be dispensed 
with, arguing that the strength of its interaction is weak. The coupled set of equations for the EYMH system 
lead to non-trivial consequences and provide a fertile ground for the study of the interplay between gravitation 
and other interactions. 

The present paper is an extension of the work by Nguyen and Wali |l2| to include electric charge and 
to study dyons coupled to gravity. The starting point is a modified EYMH system with a specific coupling 
of the Higgs field to the Einstein term in the action. In the static case, this enables us, with the help of a 
Bogomol'nyi-type [9] condition, to reduce the energy functional to a form that resembles the energy functional 
in flat space-time and derive a lower bound on the energy and hence the mass of the dyon in terms of its 
electric and magnetic charges. 

In the next section, we begin with a review of the general formalism and the field equations for the 
coupled EYMH system. In Section 3, we derive the field equations and the energy functional in the static 
case. Through the Higgs field equation, we express the Einstein term in the action in terms of metric fields, 
find a positive definite expression for the energy functional, and derive a lower bound on the energy. We also 
discuss in this section, the relation between the mass and the charge of the dyon. Section 4 is devoted to the 
derivation of the basic equations in the context of a spherically symmetric static metric and the spherically 
symmetric 't Hooft-Polyakov ansatz for the gauge and Higgs fields. In Section 5, we specialize to the Higgs 
vacuum and find monopole and dyon solutions. The final section is devoted to some concluding remarks. 



2 General Framework; Field Equations 

We begin by defining the action, 

S = J d 4 x^(£ E + £ M ), (2.1) 

where 

R-2A 9 , s 

with R, the Ricci scalar, A, the cosmological constant and the Higgs scalar field. Our metric is chosen 

to have signature (H ) with indicies fi, v, . . . running from to 3 and indicies . . . from 1 to 3, also 

q = det I 9fw\- The matter content is given by 

£m = -^V A F^ F pA + ~<r (Z? M *) - \ (* 2 - v 2 ) 2 > ( 2 - 3 ) 

where 

F^ = --[V^,V„], (2.4) 

MX 

with 

C^ = V |i -ia[A M ,-]. (2.5) 

Thus is the field strength associated with the gauge field A M , a being the strength of the gauge 
coupling. V M is the covariant derivative with the metric compatible, torsion free connection coefficients. 
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More explicitly, 
and in the component form, 



F pi , — c^A, y d u A^ zq?[A p , A y J 



(2.6) 



^ = MS - MS + a/ ahc <^, (2.7) 

where / abc are the structure constants of a gauge group Q, which for the most part for our purposes will be 
SU(2). The scalar field <fr belongs to the adjoint representation of Q. Its covariant derivative is given by 



- ia[A M ,#] 
-ia[A M ,#], 



and in component form 



V^ a = d^ a + af abc A b ^ 



(2.8) 



(2.9) 



We note that in the broken p hase of the gauge symmetry, when the Higgs field 3? assumes its vacuum ex- 
pectation value, «1> 2 = v 2 , Ce in (2.2) is the conventional Einstein-Hilbert Lagrangian. Cm in (|2.3| ) represents 
the standard Yang-Mills-Higgs Lagrangian in curved space-time. 

By varying the action S with respect to A M , 3? and g ilVl we obtain the coupled Yang-Mills, Higgs and 
Einstein equations of motion: 



1 



-.v„ 



R-2A 



8nGv 2 



\(<f> 2 -v 2 



and 



G — R^v 



R-2A 



$2 



where the energy-momentum tensor is given by 



(2.10) 
(2.11) 

(2.12) 
(2.13) 



The terms involving in ( 2.13| ) arise because of its presence in the modified Einstein-Hilbert action |L3| in 
eqn. (2.2). We further note that although the field equations are written in terms of the covariant derivative 
V M in 2? M , we can easily show that they can all be reduced to ordinary partial derivatives since 



^<t9hv = 



-gT" = o, 



(2.14) 



for torsion-free, metric compatible connection coefficients rW = . Therefore, henceforth, in our equations, 



(2.15) 



3 Static Field Equations; Bogomol'nyi Bound 



We are interested in static solutions to equations ( [2.10] ) - ( |2.12 ). Setting the time derivatives of all fields equal 
to zero, we find equation (2.10) reduces to two equations, 



DqF -'' + 



~g? ij ) 



ia [*,D *] 
ia [*,D J $] 



(3.1) 
(3.2) 
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The Higgs field equation (2.11) becomes 
1 



[D (V^D *) + D, (V=iD«#)] = - f + A I <1> 2 ! *. 



-5 ^ ' J \ 8ttGw 2 

In order to find a Bogomol'nyi-type first-order equation ]l4| to our problem, we make the ansatz |lf 

where w fc = d k f is an arbitrary time-independent function, and g = det \gij\- 
With the above ansatz, we find 



(3.3) 



(3.4) 



DiF y ' = ia [*,D J *] 



Substituting (3.5) in (|3.2|), we have 



[Do,F^] 



\/5oo 



/.9oo 



a j F y = o. 



For future reference, we note that the Yang-Mills equation (3.1) implies 

T>i (y^g~F t0 ) • * = 0. 



(3.5) 



(3.6) 



(3.7) 



Using (3.4) together with ( |3.5| ), after some algebra and the use of the Bianchi identity, the Higgs field equation 
.31) reduces to 



DnD°$ 



V9oo J V^g 

R-2A 
8ttGv 2 



\(<f> 2 ~v 2 ) $ 



In deriving the above equations, with static spherically symmetric solutions in mind, we have assumed 

goo > goi = 0, and - g > 
The static energy functional, £ , that follows from ( J2 . l| ) is given by 

£ = 



_ f°* • F 0J + ~FV ■ P„ - D 1 * D^i- (* 2 - « 2 ) 2 
8ttGv 2 2 J 2 V ; 



Next we define electric and magnetic fields, E^, and corresponding E l , B* to be given by 



E; 



E' 



/5oo 



a} 



2 2 ^-g 



and 



X = \/5oo*, 

where all the fields are functions of spacial coordinates only. 



(3.8) 
(3.9) 

(3.10) 

(3.H) 
(3.12) 

(3.13) 



With these definitions and after substituting in (3.10) for (R — 2A) /8ttGv 2 from (3.8), we obtain 

-9-, 



d 6 x<- 



E* • E,; + W ■ B, 



/3oo 



'goo 



1 V^ff / d l V9ao 

2 goo V V9oo 



z 3oo 



3oo 



3oo 



(3.14) 



In the spherically symmetric case in which we are interested, with the assumed signature for the metric, 
we note that 



.</,., 0,i t j 



and 



gu > 



(3.15) 



Consequently, each of the terms in the first parenthesis in (3.14) is positive. They are exact analogues of 
the corresponding terms in the case of flat space-time. The next two terms vanish in the Higgs vacuum, 
X 2 = goov 2 ■ Finally, the last term is a total divergence, giving rise to a finite surface term. Ignoring the terms 
that vanish when x 2 — goov 2 , we have a reduced positive-definite energy functional, £, given by 



« - \ 



d 3 



E 4 • Ei + B i • B 4 



/3oo 



-9 



3oo 



3oo 



/3oo 

\ 2 d l f 



(3.16) 



As in the case of flat space-time p6[ , we can write 



8 = 2 



d 3 x 



B 1 - cos 



E l - sin 



-9 



'goo 



9oo 

D l x ) ( Bi - cos 



D* X ) ( Ei - sin 

3 



/3oo 



/3oo 

1 



^B i .D l x+^D x-D°x+-a i 

/3oo 300 * 



vgoo 



300 



~ 9 -x 2 d*f 



(3-17) 



(3.18) 



Hence, £ , has a lower bound, 



-3 



/300 



E l ■ B.x + cost 



":B l .D !r 



-9-, 



1 



7 3oo 



^Dox • D°x + -di 

goo 2 



300 



~ 9 -x 2 d l f 



(3.19) 



reaching the lower bound when the Bogomol'nyi-type equations are satisfied, that is when 



E, - sin t 



'goo 



:T>iX = Bi - cos( 



'goo 



:D lX = 0. 



Now from our definition for B^ in ( |3.12 ), 



(3.20) 



<f~g 



B' 



1 



ijk-p 



V3oo 2 

and from the Bianchi identity, we have the conservation law, 



D, 



/3oo 



=B* I = 0. 



(3.21) 
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Consequently, 



cPx- 



-9 



/goo 



B l DiX= / d 3 xT> 



¥~9 

\/9oo 



(3.22) 



Similarly, 



/goo 



_9_ jii pviO 



(3.23) 



and 



D, 



'9oo 



E l - X ) = x-DiP i0 +E 



<T9 
\/5oo 



E' 



/5oo 



(3.24) 



on account of (3.7). 
Hence, 



-9 



'9oo 



'9oo 



(3.25) 



When the Bogomol'nyi equation ( 3 . 2C| ) is satisfied, the lower bound on the energy functional is satisfied and 
we have, 



sinfl / d i xD i 



-9, 



'9oo 



E l • x + cos 9 / d 3 xD 



Vgoo 



Bi x 



+ / d 3 x 



-9, 



1 



^D oX • D°x + -d { 

goo 2 



5oo 



~ 9 -x 2 d l f 



(3.26) 



Now, if we have finite-energy configurations with finite extension and asymptotically flat space-time, 



D^ = 



xJl 2 
<P = V , 



(3.27) 



leading to the condition that the gauge potential A M is given by 



(3.28) 



where is an arbitrary Abelian field. The field strength corresponding to the above gauge potential is 



where 



OLV A 



(3.29) 



In the static case, asymptotically, 



and 



F 40 = *-^ = -E 4 . 

v 



where we have assumed Wj = 0, i = 1,2,3. 

The above field configurations imply that outside a finite region, the non-Abelian gauge field is purely in 



the direction the direction of the unbroken U(l) electromagnetic field F^. Using ( 3.30 ), we can convert 
the divergence integrals of E 1 and B l into flux integrals of electric and magnetic fields over a surface at infinity 
and obtain 

£ = sxn6Q E + cos 6Q M + \ f d 3 xd l (—X 2 d l f) , (3.31) 

where 



2 7 \ goo 



Qe = lim / da^-id^x 2 , (3.32) 
R -^°° Js R v 



and 



Which leads to 



Qm = lim - — ^ / da m e ljk x ■ [djx, d, 
fl^oo 2av 6 J Sr 



where n, an integer, is the winding number of the mapping x '■ •S'So ~~ * & i 

Qe has no such interpretation, but its fmiteness leads to a condition on Wq. 
If we define 

. n Qe n Qm 

sin = — , „ = , cos 



>kX] ■ (3.33) 



47T 

Qm = —n, (3.34) 
a 

2 . c2 



\J Q\ + Q 2 M V Q% + Q\i 

then 



?E + Q 2 M + \ 





(^x 2 d*f) 




V 5oo / 



the equality holding when equation (3.20) is satisfied 



(3.35) 



4 Basic Equations with Spherically Symmetric Ansatz 

The considerations in the previous section are valid for any compact group Q. In this section, we shall 
specialize to SU(2). We define a spherically symmetric static metric 

ds 2 = A 2 (r)dt 2 - B 2 {r)dx 2 - r 2 (d9 2 + sin 2 6dip 2 ) , (4.1) 
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and assume the spherically symmetric 't Hooft-Polyakov ansatz 0] for the gauge and Higgs fields, 



A a = - J(r), 



„ h A b -e V (1 " W{r)) 



, nab = (-1,-1,-1), 



and 



$ a = x a vH(r). 



In spherical polar coordinates, 



(4.2) 
(4.3) 

(4.4) 



Aa _ 
A Q — 



A a n = 



J(r) 



1 - W(r) 



sm a cos ip 
sin 9 sin ip 
cos 9 
sin^ 
- cos ip 




Ai = 0, 



Aa _ 



1 - W(r) 



cos a cos ip 
cos sin ip 
— sin 9 



leading to the field strengths 
J'(r) 



pa 



17a 



sm a cos </? 
sin 9 sin (/? 

COS0 



J{r)W{r) 



a 



— sin 9? 
sin sin 99 
cos# 



pa 

rip 



W'{r) 



COS ft COS ip 

cos sin <£> 
— sin 9 



sm f 



pa 



pa _ 



J{r)W{r) 



W'(r) 



a 



W 2 (r) - 1 



COS U COS ip 

cos sin ip 
— sin# 



sm ip 

- COS </5 





sm 9 cos 95 
sin 9 sin 95 
cos6> 



sm( 



and 



sm a cos 95 
sin sin ip 
cos 



, D 9 $ a = vW(r)#(r) 



COS (/ COS ip 

cos sin ip 
— sin 9 



•#(r)ff(r) 



sin ip 
- cos y> 




sin 9, 



(4.5) 



(4.6) 



(4.7) 



where 'prime' denotes derivatives with respect to r. 

Expressed in spherical polar coordinates, our ansatz ( |3.4| ) that leads to the Bogomol'nyi-type equation 
takes the form 



D-x + [/'M-^, 



r z sm 9 

Mr) 1 F 
5(r) sinfl vr ' 

-^sinflF r9 . 



(4.8) 
(4.9) 
(4.10) 
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Substituting the spherically symmetric ansatz for x and F, we obtain two independent equations, 

B(r) (l- W 2 (r)) 
H'(r) + f(r)H(r) = — ^ 

W'{r) = -avW{r)H{r)B(r). 



(4.11) 
(4.12) 



The Yang-Mills field equations (3T) and (3J3), and the Higgs field equations are transformed into 



J"(r) + 



2 A'(r) B'(r) 
r ~ A(r) ~ B(r) 



)j'{r)-^P-J{r)W{rf = 0, 



J 2 (r)W(r) ~ f'(r) 



A'(r)\ A 2 (r) T 



-/'(r) 



A'(r)\ JT'(r) 



A(r) J H(r)B 2 {r) B 2 {r) \ B(r) r 



B'ir) 



A(r) J B 2 (r) 
f'(r) 



B 2 {r) 



W'{r) 



r"(r) 



R-2A 
8ttGv 2 



+ Xv 2 (H 2 (r) - 1) 



(4.13) 
(4.14) 

(4.15) 



We would like to remark that in the case of pure monopole in the Aq = gauge then J(r) = 0. Thus 



A'{r) 
Mr) 



if f(r) 

equation ( |4.15| ) becomes 



^4, the Yang-Mills equations (4.13) and (4.14) are automatically satisfied and the Higgs field 



B 2 (t 



B'(r) 2\ A'(r) (A'{r)\ 



B(r) r J A(r) \ A(r) 



R-2A 
8ttGv 2 



+ Xv 2 (H 2 (r) - l) . 



(4.16) 



From the metric ([O]), we calculate the components of the Einstein tensor, G M „, and the Ricci scalar, R, 
and find them to be 



and 



R 



Goo 

Gee 

G (pep 

2 I £2 



B 2 



\B - 1 



2B' 
r~B 



2 A' 
r A 



A A B r B 



2 ( A" 2 A' A' B 1 2B' 
B 2 { ~A + rl "II^rB 







r 2 ( 


1 - B 




9 

r 


( A" 


B 2 




sin 2 


9Gee, 



(4.17) 



(4.18) 



The re st o f the c omponents are zero. Likewise, we express the components of the energy-momentum tensor 
using (|il|- (l4"TQ|) , and find 



l^ H ' 2 + ^W 2 H 2 + -A 2 v 2 (H 2 if 
2 B 2 r 2 4 v ' 



(4.19) 



. J_J'2 + £_ J2 W 2 + }_ w ,2 _^_( W 2_ ^4 

2A 2 A 2 r 2 2r 4 1 > 
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A 



l H a - ^-W 2 H 2 - -B 2 v 2 (H 2 - lV 

2 r 2 4 v ' 



1 

27? 



A 2 B 2 
sin 2 GTgii. 



1 ~ 2 j' 2 + ±(w 2 -iy 



r 

B 2 



H' 



X 



V (H 2 - 1) 



Using ( |4.17| ) and ( |4.19D , it is straight forward to write Einstein field equations 



(4.20) 

(4.21) 
(4.22) 



In what follows, we eliminate from them their dependence on W'(r) and W(r) by using equations (4.11) and 
(4.12) and take linear combinations of the resulting equations. These are 



r \~A + ~B 



= 8irG 



v 2 (W_ 

IT 



1 J 2 

a 2 r 2 j\2fj2 



(4.23) 



~A 



~A~B 



B 



8nG 



1 J'' 



A 2 H 2 



v B 



v 3 Ba (H' + f'H) 



(4.24) 



A" 1 (A! B'\ A'B> 2 „ „ 

h - — ; h AS 2 = 4ttG 

A r V A B 1 A B 



1 J 



, o ( H' A 1 9 , (i? 2 - l) 



.(4.25) 



5 Higgs Vacuum; Monopole and Dyon Solutions 



We seek solutions to the set of coupled non-linear equations ( |4.11 )-(4.15) and (4.22 )-(4.25) in the Higgs 
vacuum, that is, the Higgs field frozen to a constant, H = 1. The equations simplify considerably. To this 
end, we introduce the dimcnsionless co-ordinate x and define dimensionless </, 



x = avr, 



av 



(5.1) 



We shall also work in units of the dimensionless coupling AirGv 2 set equal to unity and henceforth all the 
field variables are functions of x and 'prime' will denote derivatives with respect to x. We shall suppress their 
dependence on x for simplicity. 



From eqns (4.12) and (4.14), it follows that 



= A 1 J 2 B 
1 A A 2 ' 



and hence the Bogomol'nyi conditions (4.11) and (4.12) take the form 



A' 
AB 

B 



J 2 , B(l-y) 
A 2 x 2 ' 

-It 

2y' 



(5.2) 

(5.3) 
(5.4) 
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where we have defined y=W 2 
by 



Then the remaining Yang-Mills, Higgs and Einstein field equations ( 4.13 ), ( 4.15 ), ( 4.23 )-( 4.25 ) are given 

J " + {l-^-i) J, -i e2j » = »• < 5 - 5 > 



B 2 A 



B 2 - 1 A" 

5 2— 

x 2 A 



IB' A' (A' 
x~B + ~A It 



B 1 



B ( (A' 
2J A 2 V\-A 



J' 



(5.6) 



x \~A + ~B 



IB 2 

x 2 



A' 
AB 



= 2J 1 



'EL 

A 2 



AB + A 2 



(5.7) 



A" A' / B' 
~A ~ ~A [~B 



A' \ B 2 



= 2 



A 2 



2J' 



A 2 



1 - 



2_A 

Bl + A 2 



(5.8) 



T" 



A' 


_B>\ 


A 1 






A 


B 


A 







AS 



A 2 



J 2 A 1 



J 4 



2B —77T — h i? — r 



A 2 A 



A 4 ' 



(5.9) 



5.1 Monopole solutions 

Equations for th e m onopole are obtained by setting J = in the above set of equations. At t he out set , we 
note tha t eqn. (5^5) is automatically satisfied. By taking suitable linear combinations of (5J3), ((5.7]), (5.8) 
and (|5.9|), we find 



A (3B 2 + 1) = 0, 



(5.10) 



which implies A = for B to be real. Thus, there is no monopole solution for non-vanishing cosmological 
constant in our model. With A = 0, we are left with three independent equations, 



.4' 



1 



-A=~2 B ^-y^ S = ~o-> B = l + x-, 



It 

22/ 



A[_ 
A 



or equivalently, 



2xy 



B = 



x + y — 1 ' 



x + y — 1 

A' = l-y 
A x(x + y-l)' 



(5.11) 

(5.12) 
(5.13) 



Eq uation ( 5.12 ) is the well known Abel's differential equation of the second type. The two equations in 
(|5.13 ) are determined in terms of the solutions of Abel's equation. Abel's equation has no known analytical 
solution other than y = 0, in which case, we have an Abelian magnetic monopole with metric functions given 
by 



A 



B 2 



1 

J 2 " 



(5.14) 
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This Abelian monopole solution is also an extreme Reissner-Nordstrom black hole with mass M and charge 
Q given by 



M 



A-kGv 



Q 



M 



(5.15) 



In the general case, with non- vanishing y, we can solve equation (5.12) numerically. Results are shown in 
Figs 0,| and | 
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1.5 

x 



y(1)=0.0 
y 1 =0.2 
y(1)=0.4 
y(1)=0.6 
y(1)=0.8 

y(i =1-0 



2.5 



Figure 1: Plot of family of solutions to Abel's equation, i.e. y(x), for a range of initial conditions y(l) 
0,0.2,0.4,0.6,0.8,1.0. 



3 
2.5 

2 
1.5 

1 

0.5 


0.5 1 1.5 2 2.5 3 

x 

Figure 2: Plot of A(x), B{x) and y(x) against x for an extremal non- Abelian monopole (y(l) — 0). This 
solution is a black hole since we see the event horizon at x = 1 . 

We observe that we have a family of solutions that have exponentially vanishing non- Abelian components 
y = W 2 . They are determined by the initial condition y(l) = y^t (see Fig. |l|). These solutions are charac- 
terized by a non-Abelian core, outside of which only an Abelian component remains. The solutions to the 




12 





0.5 1 1.5 2 2.5 

x 



Figure 3: Plot of A(x), B{x) and y(x) against x for a non-Abelian monopole with the initial condition 
2/(1) = 0.5. 



metric components A and B in general have a minimum and maximum respectively in the vicinity of x = 1. 



From eqn. (5.13) the metric component A is dependent on an initial condition, however, once one fixes its 
asymptotic behavior (A — > 1 as x — > oo) it is uniquely determined. The metric component B automatically 
has the correct asymptotic behavior. We see that as the initial value of y at x = 1 approaches zero, the 
minimum of A also approaches zero at x = 1 (and the maximum B at x = 1 gets larger). The solution for 
y(V) = (see Fig. |^) represents the extremal case where an event horizon has formed. This extremal solution 
is a black hole with a non-Abelian magnetic monopole. 

For the extremal solution, the non-Abelian magnetic field is confined within the black hole horizon. Thus 
it is only natural that the metric coefficients outside the horizon are identical to that of a Reissner-Nordstrom 
black hole. In Fig. || we sec that the metric coefficient for A vanishes inside the horizon, x < 1. This is not 
a problem since A has been normalized at x — oo. An observer at infinity can not observe the interior of 
the black hole since for him an object takes an infinite amount of time to reach the horizon. If one choses to 
normalize A at the origin one would have a perfectly well defined metric inside the horizon which is infinite 
outside. This does not affect the determination of the non-Abelian magnetic field y which is independent of 
the normalization of A. This behaviour is also observed in the monopole solutions of Lue and Weinberg Q. 

In Fig. U we see a non-extremal monopole solution that is not a black hole. 
5.2 Dyonic solutions 



The set of seven coupled non-linear equations (5.3)-(5.9) involve four functions, y, J, A and B. We do not 



expect to find analytic solutions to this set of equations. Thus we are forced to look for numerical solutions. 



First, however we can analyse some features of these equations. Equation (5/7) is a quadratic equation in 
J 2 1 A 2 which has the solutions 
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Combining equations (5.3) and (5.16) for J 2 /A 2 , we obtain 



A' 
AB 



y 




AB J 



2 (A' 
l^B 2 \~A 



B> 
~B 



(5.17) 



By inspection, we see that if y = 0, the case of an Abelian dyon, then the above equation implies that 

(5.18) 

assuming A, B ^ 0. This in turn implies that the product AB is constant. 



A' B' 

A + B=°> 



Thus the Yang-Mills field equation (5.5) reduces to 

J" + 2 J' = 0, 

which has the analytic solution 

J(x) — a + c 2 /x, 



(5.19) 



(5.20) 



where C\,c% are arbitrary constants. 

Furthermore, if the charge distribution is assumed to vanish asymptotically, then c\ — 0. Together with 
the asymptotic conditions on the metric functions A, B, we can show that the cosmological constant A must 
vanish. However, inspection of equation (5.9) leads to a contradiction; which proves that there are no Abelian 
dyonic solutions to our equations (in contrast to the monopole case). 

For the general case of non- Abelian dyon solutions with y ^ 0, examination of the asymptotic behaviour 



of equation (5.9) reveals that for consistent solutions the cosmological constant must vanish. Proceeding, we 
find two coupled equations in y and J, whose solutions yield consistent solutions to all the equations. These 



are 



y"yx 3 - -y' 3 (x 2 - 2y) - -y' 2 (l - 7y + 2x 3 ) + 2y'yx 2 + 2yx 2 (x + y')[y' + y 



J 



(5.21) 
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5y 2 + 2y-l) y'(l-5y) 
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4.r 



= 0. 



The solutions to y and J yield the metric functions A and B via the following equations: 

2 



£ 
~A 

B 



li- 
ly 

2y 



i 



l-3y 



'y' 



j 



It 

2 y 



2 

x 



(5.22) 



(5.23) 
(5.24) 



Equation (5.22) is a quartic polynomial in J' / J in terms of y and y' . Since J is positive and the requirement 
that J vanish asymptotically implies that J'/J should be negative everywhere (as long as J is montonic). 
Thus we select the negative real root of the quartic polynomial. Equation (5.21) is a non-linear second order 
differential equation in y which can be readily integrated numerically using a 4th order Runga-Kutta with 
initial conditions set at x = 1. 
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Figure 4: Plot of y(x) and J(x) against x for a non-Abelian dyon with initial conditions y(l) — 1, j/(l) = — 1, 
A(l) = 2, J(l) = 1. 



Inspection of equations (5.21) and (5.22) reveals that in order for y to vanish, it can only occur at x = 1. 
One can perform a similar inspection of Abel's equation in the monopole case where it is also evident that 
y can also only vanish at x = 1. Thus, like in the monopole case one expects that the extremal black hole 
solution will be reach as the initial condition for y(x = 1) — > 0. However, unlike the monopole case there is 
an additional parameter, the initial value of y' at x = 1. 

Figures || and || contain plots of the solutions for y, J, A and B for initial conditions y(l) — 1, y'(l) = — 1, 
A(l) = 2 and J(l) = 1. The initial value of A is set in order to obtain it's correct asymptotic behavior. For 
these initial conditions the metric function A does not have a minimum and B does not have a maximum 
at finite x. However, as one adjusts the initial conditions closer to the critical value, y(l) = 0, a minima 
and maxima begin to appear. This is evident in figures ^ and |?] where the initial conditions arc y(l) = 0.5, 
y'(l) = -1.3, A(l) = 0.4, J(l) = 1. 

While inspecting the parameter space of initial conditions for y'{\) one finds that y'(l) € (—2, 0) in order 
for solutions to exist. For a value of y'(l) in this range one finds that y(l) can not be made arbitrarily close 
to zero if one requires the numerical solution to be continuously defined for all x. The smallest value of y(l) 
for a well defined solution depends on the value of y'(l). Solutions for values of y(l) closer to zero exhibit 
numerical singularities in the region x < 1. The probable explanation for this is that unlike an extremal 
monopole which only has only one horizon, a dyonic black hole naturally has two horizons. Thus, as the 
solution approaches the black hole limit the appearance of one of the horizons might occur first (meaning 
the metric coefficient A becomes small and B beco mes large at some value of x). Also, with the existence of 
two horizons it is clear that our metric (equation (4.1)), which has positive coefficients A 2 and B 2 , can not 
properly describe the region between two horizons where the space- like and time-like co-ordinates flip. Thus, 
it is not unreasonable that as our dyonic solution becomes a black hole (by tuning the initial conditions closer 
to their critical values), that the co-ordinates we used are not appropriate for all values of x. 



6 Conclusions 



We have shown in Section 3 that, with an ansatz, eqn. (3.4), we are led to Bogomol'nyi type equations and 
an energy functional that is bounded from below in the case of a dyon. In flat space-time, solutions that 
saturate the Bogomol'nyi bound are known to exist |16|. Demonstrating this in the case of a non-Abelian 
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Figure 5: Plot of A(x) and B{x) against x for a non-Abelian dyon with initial conditions y(l) — 1, y'(l) = —1, 
A(l)=2, J(l) = l. 



dyon in curved space-time is a new result. These results follow from the form of the action in eqn. (2.1), where 
the Higgs scalar is directly coupled to the Ricci scalar. In the monopole case, this unconventional coupling 
led to a first order classic Abel's equation, the solution of which yielded solutions to all the equations in the 
problem [fl2|| . In the dyon case, the problem turns out to be more complicated; nonetheless, contains similar 
simplifications in contrast to the more standard approaches in the literature. 

In Section 5, we specialize to Higgs vacuum solutions to the general equations in Section 4. For the 
sake of completeness, we begin with a review of the non-Abelian monopole solutions p2[ . We expand on 
the discussions therein. These monopole solutions depend upon initial conditions; for a specific choice the 
solution represents a non-Abelian extremal black hole. In addition we show that in our model there are no 
solutions with a non- vanishing cosmological constant in contrast to the conventional model in [^0[ . We also 
present numerical evidence for non-Abelian dyonic solutions. For a range of initial conditions the defining 
equations yield stable numerical results which are well-defined to the center of the dyon. However, as the initial 
conditions are tuned closer to critical values, the numerical results have singularities that likely correspond to 
the formation of horizons. Also as in the monopole case, physically interesting solutions require a vanishing 
cosmological constant, and further, unlike in the case of the monopole, there is no Abelian dyon solution in 
our model. 

Finally, to the best of our knowledge, a direct proof of Dirac charge quantization in the case of the EYMH 
system does not exist. However, there are compelling reasons to think, that the topological considerations 
in the case of flat space-time can be extended to curved space-time. If so, it has some extremely interesting 
consequences that have been pointed out by Ignatev, Joshi and Wali |Tij| . Among them, are the consequences 
that magnetically or electrically charged black holes obey a lower bound on their mass and the consequence 
that is not widely appreciated, namely, the spectrum of magnetically (or electrically) charged extremal black 
holes are evenly spaced in mass due to charge quantization. 
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Figure 6: Plot of y(x) and J(x) against x for a non-Abelian dyon with initial conditions y(l) = 0.5, y'(l) 
-1.3, A(l) = 0.4, J(l) = 1. 
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